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Abstract. Preliminary simulation and optimization studies of the Di- 
rectional Dark Matter Detector and the Directional Neutron Observer 
are presented. These studies show that the neutron interaction with the 
gas-target in these detectors is treated correctly by GEANT4 and that 
by lowering the pressure, the sensitivity to low-mass WIMP candidates 
is increased. The use of negative ion drift might allow us to search 
the WIMP mass region suggested by the results of the non-directional 
experiments DAMA/LIBRA, CoGeNT and CRESST-II. 
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[~ 1 Introduction 

'xT ■ 

\ Monte Carlo simulations are not only essential tools for the comparison of theory 
and experiment in physics, but also for the design and optimization of detectors. 
In these proceedings, we present preliminary simulation studies of Time Projection 
Chambers (TPCs), where the drift charge is amplified with Gas Electron Multi- 
pliers (GEMs) and detected with pixel electronics. This TPC configuration is 
described in these proceedings of S. Vahsen [T] and also in [2]. These technologies 
should allow improved gas-target detectors, where the ionization in the target gas 
^ ■ is detected with low noise, good position and time resolution, and high efficiency. 
^ \ These features allow us to measure the momentum and energy of charged particles 
and indirectly of fast neutrons by measuring charged recoils procured when they 
scatter elastically off the nuclei of the gas-target. TPCs with GEMs and pixels may 
also allow dark matter searches with improved sensitivity and background rejec- 
tion that exploit the predicted twelve-hour oscillation Weakly Interacting Massive 
Particles (WIMPs), which results from the Earth's rotation 
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Fig. 1. Flow chart showing the simulation steps expected for photons, WIMPs, nucleons, 
nuclei, leptons and mesons. 



We will present ongoing simulation studies, including first sensitivity-estimates 
for a Directional Dark Matter Detector (D 3 ) and a Directional Neutron Observer 
(DiNO) based on these technologies. These proceedings are divided into four parts: 
the simulation strategy is introduced, the preliminary simulation validation 
of GEANT4 [5] and of the WIMP cross section limit code are discussed, the design 
optimization by varying the pressure and by using either the electron drift (ED) 
or the negative ion drift (NID) is presented, and finally the preliminary results 
are shown. 

2 Simulation strategy 

The D 3 /DiNO Monte Carlo simulation toolkit is currently being assembled and 
will mostly rely on already available simulation programs: GEANT4 [5], SPJM 
[6] and GARFIELD [7]). The flow chart in Figure 1 shows the simulation steps 
expected, which will perform the following tasks: 

• An event generator will simulate the signal and background sources. E.g., 
for the WIMP case: the predicted WIMP velocity distribution for the signal 
and cosmic-ray-induced radiation and radiation emitted by detector material 
for the background. 

• An accurate geometry and materials description. 

• The interaction between the incoming particle and the target gas or/and de- 
tector materials will be modeled using either the GEANT4 [5] physic classes, 
SRIM [S], or GARFIELD [7], depending the particle type, which may lead 
to the creation of ionization along the trajectory of the incoming particle 
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Fig. 2. Left: Probability of interaction, as a function of the neutron kinetic energy, 
between the neutron and iso-butane gas-nuclei (red point and black square by isolating 
the inelastic scattering contribution) and between the neutron and the individual atomic 
nuclei (blue triangle up hydrogen, green triangle down carbon). Right: Spin dependent 
cross section limit as function of the WIMP mass for DMTPC [12] (green line) and as 
estimated by our fast simulation, assuming 100 % (dashed green line) and 50 % (red line) 
light collection efficiency. 



or/and along the recoil product of the incoming particle scattering off gas- 
nuclei. SRIM[6 and GARFIELD [7] can model theses processes fairly well. 



• Then the electrons [8] (or negative ions [9]) drift under the influence of the 
electric field. Negative ions can be formed in the case the electrons from 
the ionization can attach themselves to the gas-molecule to form negative 
ions). GARFIELD or a fast Monte Carlo simulation using the gas properties 
calculated by MAGBOLTZ [10] can simulate the drift of ionization (electrons, 
or negative ions in the case of negative ion drift), towards the GEMs with 
a constant velocity in a homogeneous electric field. The GEMs which then 
amplify the signal. In an area of high field near the GEMs, the electrons 
detach from the negative ion. Therefore a normal avalanche occurs both in 
the case of electrons and negative ions. We plan to model the GEMs with a 
parameterized simulation. 



The digitization software then simulates how the resulting avalanche-charge 
is detected by the electronic readout, in our case pixel electronics [2] ■ 
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3 Simulation validation 

3. 1 Neutron 

The probability of interaction between a neutron and the gas-target (50 cm length 
of iso-butane at 1 atm) has been simulated as function of the neutron kinetic 
energy (see Figure 1 left) with the transport code GEANT4 [5]. This probability 
has also been calculated for 1-MeV neutrons by using the Low Energy Nuclear 
Data (LEND) [IT]- There is 0.11 % per centimeter probability at 1 atm of 1-MeV 
neutron interacting with a Hydrogen atom belonging to the iso-butane-molecule. 
This analytical calculation is in good agreement with the GEANT4 [5] calculation, 
which also uses the LEND [TT] for the cross section values, and gives 0.114 % per 
centimeter. 

3.2 WIMP 

A fast Monte Carlo simulation for estimating WIMP cross section limits has been 
implemented following the instructions of Lewin et al. [3] , using SRIM [6] for the 
track length simulation and MAGBOLTZ [10] for the target-gas properties. The 
results have been compared to the DMTPC [H] published limit by applying the 
DMTPC setup parameters. A fair agreement is found, as illustrated by Figure 2 
(right). 

4 Design optimization 

The design optimization consists of finding the optimum pressure for a directional 
dark matter detector, D 3 , when there is a good trade off between the target mass 
and track length so that the directional sensitivity is maximized. The volume is 
kept fixed by considering a detector with a one square meter readout plane and a 
maximum drift length of 33.33 cm. 

Several conditions are imposed on the track to ensure that the projection of 
the track, with a length L, on the two-dimensionally segmented pixel chip readout 
plane can be exploited to extract the directionality: 

• L > 3 <j xy where a xy is the transverse diffusion 

• L > 3 x GEM holes spacing 

• energy threshold on the primary ionization energy of 1-keV that corresponds 
approximately to 40 electrons detected 

An approximate quenching factor, estimated with SRIM, was also included. 

ED and NID, when the ratio between the electric field and pressure is kept 
constant, the transverse diffusion scales with the square root of the pressure (or 
the number density). Similarly, the Townsend and attachment coefficients scale 
with the pressure. Therefore if we decrease the pressure the transverse diffusion 
will get worse and the probability of being in a multiplication regime increases 




Fig. 3. Fluorine in CF4(- top row) and sulfurine in CS2 (-bottom row): recoil (kinetic) 
energy distribution (first column), recoil track length distribution (second column), track 
length versus drift distance for fO GeV/c 2 WIMP mass (third column), angle between 
the incoming WIMP and the scatter recoil versus recoil (kinetic) energy for 10 GeV/c 2 
WIMP mass (fourth column). The directionality conditions are applied. For the recoil 
energy and track length distributions, the distributions for 100 GeV/c 2 WIMP mass are 
scaled by a factor 5. 



as well. However, the track length scales with the pressure. If we decrease the 
pressure the track length increases. 

The figure of merit (FOM) is given by the number of WIMP recoils expected 
to result in a reconstructiblc track. It has been calculated for several gases (H2, 
C 2 H 6 , C4H10, 4 He, CF 4 , 40 Ar and 132 Xe for the ED case and CS 2 for the NID 
case) and for two WIMP masses (10 GeV/c 2 and 100 GeV/c 2 ). The figure of merit 
can be written as: 



dT R Ufa dT R L{T R ,P) 

where: 

• fiA and fix are the nucleus and the nucleon reduced mass, respectively 

• p is the number of gas-molecules per cubic centimeter 

• V is the target gas volume 

• T A — F 2 (qr n )I with F 2 (qr n ) the form factor as defined in [3] and I=A 2 for 
the spin independent case (SI) 01 I k J(J+1) for the spin dependent case 
(SD) 
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Fig. 4. Figure of merit versus pressure for CF4 (left) and CS2 (right) calculated respec- 
tively for 10 GeV/c 2 (red, full line) and 100 GeV/c 2 WIMPs (black, dashed line which 
is scaled by factor 1350 for CF 4 and 450 for CS2). 



• L is the track length and Lo is the track length "threshold" derived by the 
three conditions described above. Lo is not a constant and is changing for 
each pressure 

• P is the pressure 

• Tr : is the kinetic energy of the recoil nucleus 

Example distributions (track length and recoil energy) corresponding to the 
optimum pressures found for WIMP masses of 10 GeV/c 2 and 100 GeV/c 2 , after 
the directionality conditions are applied, can be seen in Figure 3. Figure 3 also 
shows the track length versus the drift distance and the angle between the incoming 
WIMP and the recoil nucleus versus the recoil energy. For the 10 GeV/c 2 WIMP 
mass case, the most probable track length and recoil energy are: 15 mm, 13 keV 
and 0.8 mm, 5.8 keV for respectively fluorine in CF4 and sulfurine in CS2. As 
expected in the case of NID, the sulfurine recoil in CS2 have fairly short track at 
the optimum pressure. The cutoff at 0.7 mm in the track length distribution for 
CS2 corresponds to the GEM holes spacing requirement. 

Finally, Figure 4 shows the resulting figure of merit for fluorine recoils in CF4 
and sulfurine recoils in CS2. The optimum pressures to detect a 10 GeV/C 2 WIMP 
mass are 1 torr and 3 torr, respectively for CF4 and CS2. It appears that only 
CS2 among all the gases tested is competitive and can run at such low pressure. 
For example, MAGBOLTZ [TO] indicates that at 1 torr CF4, one will be in a 
multiplication regime for an electric field of few 10's V/cm, above 100's V/cm at 
1 torr one is far away the thermal limit. The gases with light nuclei can also run 
for the low pressure given by the FOM but are not competitive. For a 100 GeV/c 2 
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Fig. 5. D 3 cross section limit as a function of the WIMP mass for one recoil produced 
by a WIMP detected in three m 3 . The detector is divided into nine sub-detectors with 
a maximum drift distance of 33.33 cm for ED-CF4 and NID-CS, the SI case on the left 
and for the SD case on the right. The D 3 reach plot is compared to the non-directional 
experiments DAMA/LIBRA [13], CoGeNT 14 and CRESST-II [15] for the SI case and 
to COUPP [16] for the SD case. 



WIMP mass only CF4 and CS2 are competitive. 

For CS2-NID, a very low transverse diffusion of 50 [im/i/cm at 80 torr and an 
electric field of 1 kV/cm have been reported in [5J and is expected to run also at 
very low pressure with an electric field of few 10's V/cm according to [9]. While 
for comparison CF4 has a transverse diffusion of 150 fim/y/cm at 80 torr and an 
electric field of 2 kV/cm. 

5 Preliminary reach plots 

Figure 5 shows preliminary reach plots for WIMP search with three m 3 of CS2 or 
CF4, for three years of exposure. The detector is divided into nine sub-detectors 
with a maximum drift length of 33.33 cm. The WIMP velocity distribution used is 
Maxwellian, the dark matter density is 0.3 GeV/c 2 /cm 3 , and the escape velocity 
is 530 km/s. 

The sensitivity to low mass WIMPs is increased for low pressure as illustrated 
by the CS2 results which can be seen in Figure 5 (left). If the pressure is increased 
at the cost of directionality, the sensitivity is drastically improved as expected 
(see D 3 -ED, P = 40 torr no directionality curve in Figure 5 - left) so that after 
only a month of exposure the reach will be already below he region suggested 
by DAMA/LIBRA [13] and CoGeNT [14]. In principle a carefully designed TPC 
could decrease the pressure or even change the gas to CS2 , once a possible WIMP 
candidate is detected. Then by switching to directional mode, after three years 
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DAM A/LIBRA [13], CoGeNT [14] and CRESST-II [15] can be either excluded or 
confirmed unambiguously. 

6 Conclusion 

Preliminary simulation and optimization results of the Directional Dark Matter 
Detector and the Directional Neutron Observer have been presented. These results 
show that the neutron interaction with the gas is treated correctly by GEANT4 [5] 
and that by lowering the pressure, the sensitivity to low-mass WIMP candidates is 
increased. The use of negative ion drift might allow us to search the WIMP mass 
region suggested by results of the non-directional experiments DAMA/LIBRA |13j . 
CoGeNT rjj] and CRESST-II [15]. 

References 

[1] S.E. Vahsen, The Directional Dark Matter Detector (D 3 ) CYGNUS 2011 proceedings. 

[2] T. Kim et al, NIM A 589 (2008) 173 @ S184. 

[3] J.D. Lewin, F.P. Smith, Astroparticle Physics 6 (1996) 87-112. 

[4] S. Ahlen, N. Afshordi, J. Battat, J. Billard, N. Bozorgnia, et al., The 
case for a directional dark matter detector and the status of current ex- 
perimental efforts, Int.J.Mod.Phys. A25 (2010) 1 @ S51. larXiv:0911.0323l 
doi:10.1142/S0217751X10048172. 

[5] "GEANT4" geant4.cern.ch. 

[6] "SRIM" M. D. Z. James F. Ziegler, Jochen P. Biersack, SRIM: The Stopping and 

Range of Ions in Matter, Lulu Press Co., 2009. srim.org. 
[7] "GARFIELD" R. Veenhof, GARFIELD, recent developments, Nucl.Instrum.Meth. 

A419 (1998) 726 @ S73 0. doi:10.1016/S0168-9002(98)). 008511. 

http://garfield.web.cern.ch/garfield/ 



[8] F. Sauli, Drift and diffusion of electrons in gases: a compilation (with an introduction 

to the use of computing programs) CERN, 1984. - 127 p. 
[9] C. J. Martoff, Negative Ion TPC for WIMP AstronomySNIC Symposium, Stanford, 

California - 3-6 April 2006. 
[10] |http://consult.cern.ch/writeup/magboltz/ 



[11] National Nuclear Data Center, http://www.nndc.bnl.gov/exfor/endf00.jsp 



[12] "DMTPC" The DMTPC Detector. G. Sciolla et al., e-Print larXiv:0811. 2922 

[13] "DAMA/LIBRA" R. Bernabei, P. Belli, F. Cappella, R. Cerulli, C. Dai, et al., New 

results from DAMA/LIBRA, Eur.Phys.J. C67 (2010) 39 @ S49. larXiv:1002.1028l 

doi:10.1140/epjc/sl0052-010-1303-9. 
[14] "CoGenT" C. Aalseth, P. Barbeau, J. Colaresi, J. Collar, J. Leon, et al., Search for 

an Annual Modulation in a P-type Point Contact Germanium Dark Matter Detector, 

Phys.Rev.Lett. 107 (2011) 141301. larXiv:1106.0650l 
[15] "CRESST-II" G. Angloher, M. Bauer, I. Bavykina, A. Bento, C. Bucci, et al., Results 

from 730 kg days of the CRESST-II Dark Matter Search, larXiv: 1109.07021 
[16] "COUPP" Behnke et al., Science 319 p. 933 (2008). 



